Metal-insulator crossover in the Boson-Fermion model in infinite dimensions 
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The Boson-Fermion model, describing a mixture of tightly bound electron pairs and quasi-free 
electrons hybridized with each other via a charge exchange term, is studied in the limit of infinite 
dimensions using the Non-Crossing Approximation within the Dynamical Mean Field Theory. It 
is shown that a metal-insulator crossover, driven by strong pair fluctuations, takes place as the 
temperature is lowered. It manifests itself in the opening of a pseudogap in the electron density of 
states, accompanied by a corresponding effect in the optical and dc conductivity. 
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The discovery of high-T c superconductors (HTS) has 
led to questioning a large number of basic physical con- 
cepts. In particular, great attention has been devoted 
to the conjecture that in the normal state HTS do not 
behave as standard Fermi liquids pj. A signature for 
the possible breakdown of this concept might be the ex- 
istence of a pseudogap in the normal state which opens 
up at a temperature T* which, depending on doping, 
can be well above the superconducting transition tem- 
perature T c . The origin of this pseudogap is presently 
an issue of controversy. The fact that, upon lowering 
the temperature, this pseudogap smoothly merges into 
the superconducting gap in the entire Brillouin zone , 
strongly favors the idea that it is related to strong pair 
correlations which set in below a certain characteristic 
temperature T* and lead to condensation below T c . This 
does not mean that these pair correlations are due to 
some superconducting fluctuations but rather that they 
are a precursor to the superconducting state. They can 
exist without superconductivity and can in principle be 
stable all the way down to the lowest temperatures if the 
superconducting state is inhibited for some reason. If 
not, they will lead to condensation in some kind of su- 
perfmid state. Local probe tunneling spectroscopy inside 
and outside vortex cores at temperatures below T c repre- 
sents precisely this situation H. Such experiments rule 
out a pseudogap due to superconducting fluctuations and 
moreover confirm that the normal state pseudogap goes 
over continuously into a true superconducting gap as one 
moves across the border of the vortex cores. 

We have in the past attempted to understand the 
physics of HTS in terms of a phenomenological model 
consisting of two components: localized pairs of electrons 
(Bosons) and quasi- free electrons (Fermions), hybridized 
with each other via a charge exchange term. It is under- 
stood that the electrons are confined to the CuO-z planes 
while the localized electron pairs are of the form of bipo- 
larons and are confined to the dielectric highly polariz- 



able layers which are sandwiching the C11O2 planes [Q. 
The physics of this scenario is described by the Boson- 
Fermion model (BFM) 
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with So = D — fi and E — Ab — 2/i. Here c^J denote 
fcrmionic operators for electrons with spin a at some 
effective sites i (involving molecular units rather than 
individual atoms) and ft are hard-core bosonic opera- 
tors describing tightly bound electron pairs. The bare 
hopping integral for the electrons is given by t, the bare 
electronic half bandwidth by D, the Boson energy level 
by As and the Boson-Fermion pair-exchange coupling 
constant by g. The chemical potential fj, is common to 
Fermions and Bosons (up to a factor 2 for the Bosons) in 
order to guarantee charge conservation. 

The underlying physics of this model is that the charge 
exchange term between Bosons and free electrons induces 
local pairing amongst the latter and thus leads to a deple- 
tion of the density of states near the Fermi level which 
reflects itself in the opening of a pseudogap. Our pre- 
vious studies of this model on the basis of lowest order 
self-consistent conserving approximations bore this out 
Q and showed how the appearance of such a pseudo- 
gap manifests itself in the thermodynamic, transport and 
magnetic properties of the system || . These studies in- 
dicated that T* was intimately related to a precursor of 
a condensation of the intrinsically localized Bosons due 
to some highly non-linear feedback effect from the itin- 
erant electrons. This finding was verified independently 
by a method summing up the most diverging diagrams 
[7|. Our follow-up study of this model in the atomic limit 
[ B I (which excluded from the outset any global supercon- 
ducting state) gave the first indications that the open- 
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ing of the pseudogap was related to a metal-insulator 
crossover rather than to superconducting pair fluctua- 
tions. Both effects, of course, will determine the temper- 
ature dependence of the pseudogap, but with different 
importance depending on the temperature regime. We 
shall present in this Letter a dynamical mean field study 
of the BFM in infinite dimensions in which full use is 
made of those atomic limit results and which permits 
to incorporate the itinerancy of the electrons in a con- 
trolled fashion. We shall below sketch the main lines of 
this algorithm and then concentrate on the normal state 
properties of this model. 

In infinite dimensions a many body system can be 
treated as a purely local system coupled to a "medium" 
described by a Weiss field, with which it interacts in the 
form of an impurity Anderson problem In such a 
scheme the effective Hamiltonian for the BFM is 

H =^eo4e CT + E Q b f b + g [ c\c\b + tfc^ } 
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where and denote the original Fermion and Bo- 
son operators at the impurity site, and a denote the 
auxiliary Fermionic variables of the Weiss field. For 
these latter the energy spectrum w k and their coupling 
Vk to the impurity electrons have to be determined self- 
consistently. We see that the first line of Eq.(2) coincides 
with the local part of the Hamiltonian (1). Setting H = 
Hq + Hi , with Hj denoting the interaction term between 
the electrons and the auxiliary Fermions, we have that 
Ho can be written in the form Hq = Y] -, \n) E n (n| + 
H med , where H med = J^k a w k d{ a d kiCr and 



The interaction Hamiltonian then becomes 
Hi = ^2 Vk t F ™n d l a X mn + F° m X mn d k ^ } . (5) 
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We shall now solve the d — > oo BFM within the so-called 
Non-Crossing Approximation (NCA) in Hi. A similar 
approach was recently used in the Dynamical Mean Field 
Theory (DMFT) context || by Lombardo et al. for a 
multiband Hubbard model for perovskites |lCj and by 
Schork and Blawid for the Anderson lattice model with 
correlated conduction electrons [ pT| . For this purpose we 
introduce the resolvents R mn (z), where m can be differ- 
ent from n, defined through the Dyson equation 

Rmn{z) = S mn R° m {z) + R°m{z) ^ml(z) Rln(z) , (6) 



together with the associated spectral functions given by 
4 mn (w) = — 21mR mn {u> + id). The free diagonal resol- 
vents in Eq.(^) are R^z) = (z — E m )~ l and the self- 
energies S m /(z) are determined self-consistently through 
the equations 
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with the use of the Weiss field spectral function 
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are the eigenstates and the eigenvalues of the Hamilto- 
nian (1) in the atomic limit. Here the arrows denote the 
presence of an electron with spin up or down and the dot 
indicates the presence of a Boson on the effective impu- 
rity site. In Eqs.(|^) we have defined u 2 = 1/2 — (2£o — 
£o)/(4 7 ), v 2 = 1 - u 2 and 7 = ^/(e - E Q /2) 2 + g 2 . 

Introducing the Hubbard operators X mn = \m)(n\ we 
can rewrite the original Fermion operators as 



1 ran ^mn 



r nm ^-ran 



(4) 



where 

F 1 - 

F i . 
r 13 " 



^78 



1, fJ 



^68 



34 
1) F 24 



r 56 



U, F. 



r 46 



-u, 4 



25 



Fi 7 



Within the NCA the local Fermion Green's function 

Gf{z) = ^2^2 FmnFn'rn'{( X mn ] X m ' n ')) z (9) 



is expressed in terms of the resolvents according to the 
relations 
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where 
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and ^ = — 1 for Fermions and +1 for Bosons. The above 
set of self-consistent equations is evaluated for a Bethe 
lattice with a semi-circular density of states D(e) of width 
2D = At. In this case the DMFT self-consistency condi- 
tion takes the simple form A(w) = t 2 Af(lu) M, where 
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Ap{ui) = —2 ImGF(u+ii5) is the fermionic spectral func- 
tion. We should stress that our approach represents 
a generalization of the NCA (originally devised for the 
single- impurity Anderson model Jt2|]), suited for models 
with a total number of particles per site greater than 
one. The original NCA formulation is recovered when all 
the non-diagonal resolvents are neglected, as was recently 
done for the multiband correlated electron models stud- 
ied in Refs. [ p~0|JTT| . In general, the iterative procedure 
determining those resolvents will lead to an increasing 
number of such non-diagonal elements, the relevance of 
which will depend on the model under consideration. For 
the problem studied here, the only non-vanishing non- 
diagonal element is -R45 (= -R54). 

The results presented below have been obtained in the 
symmetric limit of the model, i.e., for £0 = Eq = (from 
now on all energies are in units of 2D). In this case the 
chemical potential is fixed in the middle of the electronic 
band, giving n F = J2a( c l c <?) = 1 and n B = {^b) = \. 
We furthermore choose g = 0.2 such that the character- 
istic temperature of the opening of the pseudogap lies in 
the regime of a few hundred degrees K. From the behav- 
ior of the density of states Ap(oj), plotted in Fig. 1 for 
various temperatures T, we can see that the pseudogap 
opens up below a characteristic temperature T* ~ 0.15 
and gradually deepens as T is reduced. This value of T* 
gives a reasonable estimate of the onset of the pseudogap 
regime, provided that the Fermions are assumed to be 
quasi-particles with a strongly renormalized bandwidth. 
We also notice that temperatures lower than those re- 
ported in the figure are not considered here, since they 
do not lead to reliable convergency rates. This is proba- 
bly a sign of the failure of the NCA at low temperatures. 

The behavior of the pseudogap is best studied in 
terms of the fermionic self-energy or its imag- 

inary part T(lu) = — 2ImEi?(z = u> + iS). (z) 
is introduced by expressing the local Green function 
in the form G F (z) = [z - s - E w (z) - EfO)]" 1 , 
where the Weiss self-energy Y, w (z) is determined from 
A(oj) = — 2ImS M '(aj + iS) or, alternatively, using the 
self-consistency condition T, w (z) = t 2 Gp(z). For an or- 
dinary Fermi Liquid, one expects that T(lu) ~ lj 2 such 
that the lifetime of the quasi-particles at the Fermi energy 
T F = l/r(0) becomes infinite. In Fig. 2a, we plot Y{uS) 
for the same parameters as in Fig. 1 and observe that the 
pseudogap is linked to a strong resonant scattering of the 
quasi-particles with the Bosons at the Fermi energy. As 
one lowers the temperature, the lifetime at the Fermi en- 
ergy tends to zero and a metal-insulator crossover takes 
place. The form of clearly shows that our system 
does not correspond to a Fermi liquid, as we can see from 
an inspection of the poles of the lattice Green's function 
Gp(z,lu) = [z-£i;-Sf(z)] _1 (ek denoting the bare elec- 
tronic dispersion). At high temperatures, T ^> T*, we 
see from Fig. 2b that the straight line representing the 
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FIG. 1. Fermion DOS for several temperatures T (all ener- 
gies in units of 2D), the deepest pseudogap corresponding to 
the lowest T. The inset shows the DOS at the Fermi energy 
as a function of T. 



inverse of the free-particle Green's function at A; = kp 
(with Ek F = 0) cuts Ke'Sp(uj) only at the Fermi energy 
uj = 0. Considering, however, that T(u>) develops a min- 
imum at lu = which tends to a constant as T is in- 
creased, we still cannot attribute this feature to a normal 
Fermi liquid behavior. Upon lowering the temperature, 
we observe an inversion of slope which results in a cut 
occurring at three distinct frequencies, with the solution 
at to — being now accompanied by a large imaginary 
part of T*p. The pseudogap is thus a consequence of a 
spectral weight transfer from the Fermi energy away to 
the wings in the density of states associated with bond- 
ing and anti-bonding two-particle states |^] . This is very 
similar to what has been obtained for the 2d Hubbard 
model at half filling p3| , where a pseudogap develops 
due to antiferromagnetic correlations. The results pre- 
sented here are also in line with the precursor pairing 
correlation ideas developed by Randcria JT3] . 
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FIG. 2. Frequency dependence of the imaginary (a) and the 
real (b) part of the Fermion self-energy for several tempera- 
tures T (same units as in Fig. 1). The lowest T corresponds 
to the highest value of ImEp(0) and the steepest slope of 
Re E F (0) , respectively. 
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The present study does not contain any superconduct- 
ing fluctuations and hence excludes the possibility of a 
transition to a superconducting state. The opening of the 
pseudogap is hence entirely related to a crossover from 
a metallic-like state at high temperatures to an insulat- 
ing one at lower temperatures. We now show that this 
behavior is also manifest in the real part of the optical 
conductivity which in the limit d — > oo takes the form 0] 



a normal state pseudogap due to uncorrelated pair fluctu- 
ations which, upon approaching T c from above, will tend 
to correlate and lead to a true superconducting gap. 
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Here the sum over momenta in the spectral function 
Ap(ek,u}) of the lattice Green's function has been ex- 
pressed as an energy integration over D[e). From the 
results presented in Fig. 3, we notice that the opening 
and gradual deepening of the pseudogap, as the temper- 
ature is reduced, is associated with a change-over in the 
dc resistivity <r(0) _1 from a linear in T decrease to an 
upturn (see the inset). We stress that the linear behav- 
ior does not arise from the temperature dependence of 
the single-particle lifetime r(0) _1 , which for T>T* be- 
comes independent on T, but is exclusively due to an in- 
trinsic non- Fermi liquid behavior even above T* . This is 
also manifest in the fact that the Fermion DOS never re- 
covers the unperturbed semicircular shape, however high 
T is. As far as the frequency dependence of <j{w) is con- 
cerned, we observe below T* a well defined isosbestic 
point around which the spectral weight is shifted from 
low to higher frequencies and at which the various curves 
for the optical conductivity, corresponding to different T, 
cross. This behavior, recently confirmed by optical mea- 
surements p5[ , is in line with our previous results || on 
the optical conductivity in the Id and 2d BFM. 

In conclusion, we have reported in this Letter on a 
metal-insulator crossover described by the BFM for HTS, 
driven by strong pair fluctuations and leading up to the 
opening of a pseudogap in the density of states. Our 
calculations are based on a generalization of the non- 
crossing approximation, within a dynamical mean field 
study in the d — *■ oo limit. The results obtained are simi- 
lar to those for strongly correlated systems were the pseu- 
dogap is due to antiferromagnetic correlations. In order 
to decide which of the mechanisms is relevant for HTS, 
a detailed knowledge is required of the k dependence of 
this pseudogap in the Brillouin zone, together with the 
temperature behavior of the incoherent part of the single- 
particle spectrum. A study for the latter, on the basis of 
a generalized BFM in the atomic limit, associates it with 
phonon shake-off effects, when the Bosons are considered 
as bipolarons Whatever the mechanism of pairing 

might be, the present study of the BFM clearly suggests 
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FIG. 3. Fermion optical conductivity for several tempera- 
tures T, the lowest T corresponding to the highest value of 
a(oj) at high frequencies (same units as in Fig. 1). The inset 
illustrates the dc resistivity (cr(0)/7r) _1 as a function of T. 
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